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The temperature-driven reorientation transition which, up 
to now, has been studied by use of Heisenb erg- type models 
only, is investigated within an itinerant-electron model. We 
consider the Hubbard model for a thin fcc(lOO) film together 
with the dipole interaction and a layer-dependent anisotropy 
field. The isotropic part of the model is treated by use of a gen- 
eralization of the spectral-density approach to the film geom- 
etry. The magnetic properties of the film are investigated as a 
function of temperature and film thickness and are analyzed 
in detail with help of the spin- and layer-dependent quasipar- 
ticle density of states. By calculating the temperature depen- 
dence of the second-order anisotropy constants we find that 
both types of reorientation transitions, from out-of-plane to 
in- plane ("Fe-type") and from in-plane to out-of-plane ("Ni- 
type") magnetization are possible within our model. In the 
latter case the inclusion of a positive volume anisotropy is 
vital. The reorientation transition is mediated by a strong 
reduction of the surface magnetization with respect to the in- 
ner layers as a function of temperature and is found to depend 
significantly on the total band occupation. 
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I. INTRODUCTION 

The large variety of novel and interesting phenomena 
of thin-film magnetism results very much from the fact 
that the magnetic anisotropy, which determines the easy 

' axis of magnetization, can be one or two orders of mag.- 

■ nitude larger than in the corresponding bulk systemsEJ. 
The reorientation transition (RT) of the direction of 
magnetization in thin ferromagnetic films describes the 
change of the easy axis by variation of the film thick- 

' ness or temperature and has been addfily studied both 
experimentallyoEJ and theoreticallyll3~E3. 

An instructive phenomenological picture for the un- 
derstanding of the RT is obtained by expanding the free 

| energy F of the system in powers of cos 9m > where 9 m 
is the angle between the direction of magnetization and 
the surface normal. Neglecting azimuthal anisotropy and 
exploiting time inversion symmetry yields: 
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The anisotropy coefficients of second (if 2) and fourth 
(K4) order depend on the thickness d of the film as well 
as on the temperature T. 

Away from the transition point usually K 2 3> K4 
holds, and, therefore, the direction of magnetization is 
determined by the sign of if 2 (if 2 > 0: out-of-plane mag- 
netization; if 2 < 0: in-plane magnetization) On this ba- 
sis the concept of anisotropy flowE3ii3 immediately tells 



us that the RT is caused by a sign change of if 2 while 
the sign of K4 mainly determines whether the transition 
is continuous (if 4 < 0) or step-like (if 4 > 0). In the case 
of a continuous transition K4 also gives the width of the 
transition region. 

From the microscopic point of view we know that 
the magnetic anisotropy is exclusively caused by two ef- 
fects, the dipole interaction between the magnetic mo- 
ments in the sample and the spin-orbit coupling: if; = 
ifj.so + if i. dip- While the dipole interaction always favors 
in-plane magnetization (if2,dip < 0) due to minimiza- 
tion of stray fields, the spin-orbit interaction can lead to 
both, in-plane and out-of-plane magnetization depending 
sensitively on the electronic structure of the underlying 
sample. The spin-orbit anisotropy is caused by the bro- 
ken symmetryo at the film surface and the substrate-film 
interface as well as by possible strainEai 2 -! in the volume of 
the film. It is worth to stress that a strong positive spin- 
orbit induced anisotropy alone opens up the possibility 
of an out-of-plane magnetized thin film. The RT must 
be seen as a competition between spin-orbit and dipole 
anisotropy. 

In many thin-film systems both thickness- and 
temperature-driven RTs are observed. Although it 
is clear by inspection of the corresponding phase 
diagramsoo that both types of transitions are closely 
related to each other, different theoretical concepts are 
needed to explain their physical origin. 

The thickness-driven RT is rather well understood in 
terms of a phenomenological separation of the spin-orbit 
induced anisotropy constant ifj jSO into a surface term 
iff so and a volume contribution Kj, so by the ansatz 
ifj,so = if^ S o + 2iff so /d. Experimentally, this seDara- 
tion seems to provide a rather consistent pictureQclliillij 
despite the fact that in some samples additional struc- 
tural transitions are presentBEJ which clearly restrict 
its validity. On the theoretical side, basically two dif- 
ferent schemes for the calculation of T = magnetic 
anisotropy constants Ijayeieen developed, semi-empirical 
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tight-binding theorieai3~ 
total-energy calculations! 2 ^ 
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,id spin-polarized ab initio 
In both approaches the 
spin-orbit coupling is introduced either self-consistently 
or as a final perturbation. However, these investigations 
still remain to be a delicate problem because of the very 
small energy differences involved. 

Neglecting the large variety of different samples, sub- 
strates, growth conditions, etc. it is useful for the un- 
derstanding of the RT to concentrate on two somewhat 
idealized prototype systems both showing a thickness- as 
well as a temperature-driven RT. 



The "Fe-type" systemJHl are characterized by a large 
positive surface anisotropy constant K% together with a 
negative volume anisotropy K\ due to dipole interac- 
tion. This leads to out-of-plane magnetization for very 
thin films. For increasing film thickness the magnetiza- 
tion switches to an in-plane direction because the volume 
contribution becomes dominatingcTLl. As a function of in- 
creasing temperature a RT from out-of-plane iron-plane 
magnetization is found for certain thicknessesoa 

In the "Ni-type" systemsliirtL3, the situation is differ- 
ent. Here the volumaxpntribution K~ is positive due to 
fct lattice distortion£3e3, thereby favoring out-of-plane 
magnetization, while the surface term is negative. 
For very thin films the surface contribution dominates 
leading to in- plane magnetization. At a critical thickness, 
however, the positive volume anisotropy fojxeiLthe system 
to be magnetized in out-of-plane directior£il"E3, until at a 
second critical thickness the magnetization switches to an 
in-plane position again caused by structural relaxation ef- 
fects. Here a so-called anomalous temperature-driven RT 
from in-plane to out-otplane magnetization was found re- 
cently by Farle et al.E-Kj. 

In this article we will focus on the temperature-driven 
RT which cannot be understood by means of the separa- 
tion into surface and volume contribution alone. Here the 
coefficients K~ and KY need to be determined for each 
temperature separately. Experimentally, this has been 
done in grea± detail for the second-order anisotropy of 
Ni/Cu(100)H3. The results clearly confirm the existence 
and position of the RT, but, on the other hand, do not 
lead to any microscopic understanding of its origin. 

To obtain more information on the temperature-driven 
RT theoretical investigations on simplified model sys- 
tems have proven to be fruitful. Despite the fact that 
in the underlying transition-metal samples the sponta- 
neous magnetization is caused by the itinerant, strongly 
correlated 3d-electrons, up to now Hei senhe rg-type mod- 
els have been considered exclusivelyE^Lj. The mag- 
netic anisotropy has been taken into account by incorpo- 
rating the dipole interaction and an uniaxial single-ion 
anisotropy to model the spin-orbit-induced anisotropy. 
Using appropriate T — second-order anisotropy con- 
stants as input parameters, both types of RTs have been 
observed within th e fram ework of a self-consistent mean 
field approximatior£3E3 as well a s b y first-order pertur- 
bation theory for the free energjapL A continuous RT 
has been found for d > 3 layersHE2l taking place over a 
rather small temperature range. Step-like transitions oc- 
cur as an exception for special parameter constellations 
only. The RT is attributed to the strong reduction of the 
surface-layer magnetization relative to the inner layers for 
increasing temperature leading to a diminishing influence 
of the surface anisotropy. 

Since the itinerant nature of the magnetic moments is 
ignored completely in these calculations, it is interesting 
to compare these results with calculations done within 
itinerant-electron systems. The present work employs a 



similar concept biit in the framework of the single-band 
Hubbard modelEU which we believe is a more reason- 
able starting point for the description of temperature- 
dependent electronic structure of thin transition-metal 
films. 

The paper is organized in the following way: In the 
next section we define our model Hamilton operator. In 
Sec. Ill we will focus on the derivation of the free en- 



ergy and the second-order anisotropy constants by use 
of a perturbational approach. The isotropic part of the 
Hamilton operator is treated in Sec. [fy| Here we present 
a generalization of a self-consistent spectral-density ap- 
proach to the film geometry. In Sec. ^ we will show and 
analyze the results of the numerical evaluations and dis- 
cuss the possibility of a temperature-driven RT within 
our model system. We will end with a short conclusion 
in Sec. 0. 

II. DEFINITION OF THE HAMILTON 
OPERATOR 

The description of the film geometry requires some 
care. Each lattice vector of the film is decomposed into 
two parts 



(2) 



Ri denotes a lattice vector of the underlying two- 
dimensional Bravais lattice with N sites. To each lattice 
point a (i-atom basis r a (a — 1 , . . . , d) is associated re- 
ferring to the d layers of the film. The same labeling, of 
course, also applies for all other quantities related to the 
film geometry. Within each layer we assume translational 
invariance. Then a Fourier transformation with respect 
to the two-dimensional Bravais lattice can be applied. 

The considered model Hamiltonian consists of three 
parts: 



H — -ff o + -ffdip + H so . 
Hq denotes the single-band Hubbard model 



(3) 
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where Ci aa {c\ a<y ) stands for the annihilation (creation) 
operator of an electron with spin a at the lattice site 
R-ia, ni a = c\ aa Ci acr is the number operator and T^ 3 is 
the hopping- matrix element between the lattice sites Ri Q 
and Rj/3. The hopping- matrix element between nearest 
neighbor sites is set to —t. U denotes the on-site Coulomb 
matrix element, and \i is the chemical potential. 

The second term f/dip describes the dipole interaction 
between the magnetic moments on different lattice sites: 



-ffdip — — 



1 



i,j,a,/3 (~^) 3 



(5) 



Here r"? is the distance between H ia and Rj/3, and 



the unit vector is given by u"' 3 



c^o = A i oMs z «c/(47ra 3 ) denotes the strength of the dipole 
interaction, a the lattice constant and z MC the num- 
ber of atoms in the corresponding bulk cubic unit cell. 

cr TT iCi aT i are the spin operators con- 
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structed by the Pauli spin matrices <x. The expectation 
value of the spin operator yields the (dimensionless) mag- 
netization vector 



m Q — {(Ti a ). 



(6) 



The magnetization only depends on the layer index a 
because of the assumed translational invariance within 
the layers. 

In our approach the uniaxial anisotropy caused by the 
spin-orbit coupling is taken into account phenomenologi- 
cally by an effective layer-dependent anisotropy field cou- 
pled to the spin operator: 



(7) 



The effective field Ba°' is chosen to be parallel to the 
film normal n. To ensure the right symmetry we set: 



=/3 Q (m a n)-n. 



(8) 



This corresponds to a meaii-field treatment of the 
spin-orbit-induced anisotropyEx The strengths of the 
anisotropy fields j3 a enter as additional parameters and 
have to be fixed later. 



III. CALCULATION OF THE MAGNETIC 
ANISOTROPY CONSTANTS 

The direction of the magnetization is determined by the 
minimal free energy F. The anisotropic contributions H so 
and -f/dip to the Hamilton operator can be considered 
as a small perturbation to the isotropic Hubbard model 
((3 a , ljq <C i, U). Then_wa can apply a thermodynamic 
perturbation expansionEJEj of the free energy F up to 
linear order with respect to j3 a , ujq: 

F{T) = F (T) + (H so ) /(Nd) + {H dip ) /{Nd). (9) 

Here (. . .)o denotes the expectation value taken within 
the unperturbated Hubbard model Ho. On the same foot- 
ing we use a mean-field decoupling for the two-particle ex- 
pectation values contained in (Hdip)o- Because F is cal- 
culated to linear order in the anisotropy contributions, 
only the lowest, i.e. second-order anisotropy constants 
are considered, and a possible canted phase is neglected. 
The ratio K4/K2 and thus the width of the transition 
region have been found to be very small_fpr reasonable 
strengths of the anisotropy contributionsE§E3. 

Within our approach it is, therefore, sufficient to con- 
sider the free-energy difference K2{T) between in-plane 
and out-of-plane magnetization: 



K 2 (T) = F(T,9 



-) 



F(T,6 M = 0). 



(10) 



K2 < and K2 > indicate in-plane and out-of-plane 
magnetization respectively. Hence, the reorientation tem- 
perature Tr is given by the condition ^(Tfl) = 0. Eval- 
uation of Kq,(T) yields: 



K 2 (T) = K 2 , s 



^2,dip 

2 3o;o 1 x - 



a/3 



m a m/3$\ a _i3\. (11) 



The constants $ Q contain the effective dipole interaction 
between the layers and can be calculated separately: 



*|a-/3| = — - 
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(12) 

@iaj/3 ^ S ^ ne an gl e between u^ 8 and the direction of the 
magnetization. The <E> Q only depend on the film geome- 
try. For thick films -f^.dip reduces to its continuum value 
^HoM 2 where M is the magnetization per atom. To cal- 
culate K2 (T) the temperature- and layer-dependent mag- 
netizations m„ of the Hubbard film are needed. 



IV. SPECTRAL-DENSITY APPROACH TO THE 
HUBBARD FILM 

In this section we will focus on the evaluation of 
the Hubbard film Hq. Ferromagnetism in the Hubbard 
model is surely a strong-coupling phenomenon. The ex- 
istence of ferromagnetic solutions was recently proven 
in the limit of infinite dimensions by quantum Monte- 
Carlo calculationsl3t3. Ferromagnetism is favored by a 
strongly asymmetric Bloch density of states (BDOS) and 
by a singularity at the upper band edge as it is found, e. 
g., for the fee lattice. 

The Hubbard model constitutes a highly non-trivial 
many-body problem even for a periodic infinitely ex- 
tended lattice. Even more complications are introduced 
when the reduction of translational symmetry has to be 
taken into account additionally. One of the easiest possi- 
ble approximations to treat a Hubbard film is a Hartree- 
Fock decoupling, which has been applied previousljCj'tll. 
Hartree-Fock theory, however, is necessarily restricted 
to the weak-coupling regime and is known to overesti- 
mate the possibility of ferromagnetic order drastically. 
Neglecting electron-correlation effects altogether leads to 
qualitatively wrong results especially for intermediate 
and strong Coulomb interaction Ue3. Furthermore, we 
did not find a realistic temperature-driven RT within this 
approach. 

Thus we require an approximation scheme which is 
clearly beyond the Hartree-Fock solution and takes into 
account electron correlations more reasonably. On the 



other hand, it must be simple enough to allow for an ex- 
tended study of magnetic phase transitions in thin films. 
For thi^pnrpose we apply the spectral-density approach 
(SDA)c3e3 which is motivated by the rigorous analysis 
of Harris and Langetj in the limit of strong Coulomb 
interaction (U/t 3> 1). The SDA can be intemmted 
as an extension of their tj [/-perturbation theoryc3c3 in 
a natural way to intermediate coupling strengths and 
finite temperatures and has been discussed in detail 
for various .three- di m en si on a MM-pI as well as infinitc- 
dimensionalE3'E3 lattices. At least qualitatively, it leads to 
rather convincing results concerning the magnetic prop- 
erties of the Hubbard model. A similar approach has 
been applied to a multiband Hubbard model with sur- 
prisingly accurate results for the magnetic key-jquantities 
of the prototype band ferromagnets Fe, Co, NiE3Ed. Re- 
cently, a generalization of the SDA has been proposed 
to deal with the, modifications due to reduced transla- 
tional symmetrytHrEJ. In the following we give only a 
brief derivation of the SDA solution and refer, tha mader 
to previous papers for a detailed discussionci2e5E3& 
The basic quantity to be calculated is the retarded 



single-electron Green function Gff a (E) = ( (c 



From Gf J ^ a .(E) we obtain all relevant information on the 
system. Its diagonal elements, for example, determine the 
spin- and layer-dependent quasiparticle density of states 
(QDOS) p aa (E) = -ilmGg°(£ - fi). The equation of 
motion for the single-electron Green function reads: 



' "'<ri C 'jl3<j)) E - 



E^ + ^S 7 



rpa~t 



KliE)]GtiE) 



Uf. (13) 



Here we have introduced the electronic self-energy 
Sy^. (E) which incorporates all effects of electron correla- 
tions. We adopt the local approximation for the self- 
energy which has been tested recently for the case of 
reduced translational symmetr If we assume trans- 
lational invariance within each layer of the film we have 



^ t]a (E) = 5fPY%(E). After Fourier transformation with 
respect to the two-dimensional Bravais lattice the equa- 
tion of motion ( |l3"| ) is formally solved by matrix inversion. 

The decisive step is to find a reasonable approximation 
for the self-energy (E) . Guided by the exactly solvable 
atomic limit of vanishing hopping (t = 0) and by the 
findings of Harris and Lange in the strong-coupling limit 
(U/t 3> 1), a one-pole ansatz for the self-energy H"(E) 
can be motivatedEl The free parameter of this ansatz are 
fixed by exploiting the equality between two alternative 
but exact representations for the moments of the layer- 
dependent quasiparticle density of states: 



1 



dE(E - p) m p aa (E) 



(ih 



dt> 



«r(*),4«r(*') 



(14) 



Here, [...,...] + denotes the anticommutator. It can 
shown by comparing various approximation schemes^ 
that an inclusion of the first four moments of the QDOS 



(rn = — 3) is vital for a proper description of ferro- 
magnetism in the Hubbard model. Further, the inclusion 
of the first four moments represents a necessary condi- 
tion to be consistent with the i/[/-perturbation theoryca. 
The Hartree-Fock approximation recovers the first two 
moments, (m = 0, 1) only, while the so-called Hubbard-I 
solutionEH reproduces the third moment (m = 2) as well, 
but is well known to be hardly able to describe ferro- 
magnetism. Taking into account the first four moments 
to fix the free parameters of our ansatz we end up with 
the SDA solution which is characterized by the following 
self- energy: 



K(E) = Un a - a 



E + fi — B a _ a 



E + (1- B a - a — U(l — Uct-a) 



(15) 



The self-energy depends on the spin-dependent occupa- 
tion numbers n aa = (c\ aa Ci aa ) as well as on the so-called 
bandshift B a(T which consists of higher correlation func- 
tions: 



B n 



rpcta _i_ 
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E ( c L<x c j73a(2n;a- 



3,0 



(16) 



A possible spin dependence of B aa opens up the way 
to ferromagnetic solutions00. Ferromagnetic order is 
indicated by a spin-asymmetry in the occupation num- 
bers n a i ^ n a ^, and the layer-dependent magnetization 
is given by m a = n Q | — n Q j,. 

The band occupations n aa are given by 

+oo 

n aa = J dEf_(E) Paa (E), (17) 

— oo 

where f—(E) is the Fermi function. The mean band oc- 
cupation n is defined as n — i ^2 acr n acr . Although B a _ ff 
consistsoLhigher correlation functions it can by expressed 
exactlyH@ via p aa (E) and 

B a<J =T% a + -i -i / dEf_(E)x 

n aa (l - n aa ) h J 



-Y%{E - M ) - lj [E - K(E - M) - TS a ] Paa (E). ( 

Equations @, @, @ and @ build a closed set of 
equations which can be solved self-consistently. 

Surely a major short-coming of the SDA is the fact that 
quasiparticle damping is neglected completely. Recejitly 
a modified alloy analogy (MAA) has been proposecE3£3 
which is also based on the exact results of the t/U- 
perturbation theory but is capable of describing quasipar- 
i=tticle damping effects as well. For bulk systems it has been 
found that the magnetic region in the phase diagram is 
significantly reduced by inclusion of damping effects. On 
the other hand, the qualitative behavior of the magnetic 
solutions is very similar to the SDA. An application of 
the MAA to thin film systems is in preparation!!! 



V. RESULTS AND DISCUSSION 

The numerical evaluations have been done for an 
fcc(lOO) film geometry. In this configuration each lattice 
site has four nearest neighbors within the same layer and 
four nearest neighbors in each of the respective adjacent 
layers. We consider uniform hopping T?j = —t between 
nearest neighbor sites R< a , R^g only. Energy and tem- 
perature units are chosen such that t = 1. The on-site 
hopping integral is set to T-$ a = 0. Further, we keep 
the on-site Coulomb interaction fixed at U — 48 which is 
three times the band width of the three-dimensional fcc- 
lattice and clearly refers to the strong-coupling regime. 



/ dE(E - T°l a ) n pS\E) of the BDOS yields that 
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FIG. 1. Layer-dependent Bloch density of 
states Pa (E) = p a a{E)\u=o of a two, three and five layer 
fcc(100) film, a denotes the different layers of the film where 
a — 1 corresponds to the surface- layer. The nearest-neighbor 
hopping is set to t = 1. 

Let us consider the isotropic Hubbard- film first. There 
are three model parameters left to vary, the temperature 
T, the thickness d and the band occupation (0 < n < 2). 
Except for the last part of the discussion we will keep 
the band occupation fixed at the value n — 1.4 and focus 
exclusively on the temperature and thickness dependence 
of the magnetic properties. 

In Fig. |l| the layer-dependent density of states of the 
non-interacting system pa(E) = p a a(E)\u = Q (= "Bloch 
density of states" (BDOS)) is plotted for a two, three 
and five layer film with an fcc(100) geometry. The 
BDOS is strongly asymmetric and shows a considerable 
layer-dependence for d > 3. Considering the moments 



(2) 

the variance A„ 



(3) 

as well as the skewness A a are re- 
duced at the surface layer compared to the inner layers 
due to the reduced coordination number at the surface 
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FIG. 2. Layer magnetization m a and mean magnetization 
m — 2 X/q ma as a function of temperature for different num- 
ber of layers d = 2,3,4,5. Further parameters: n = 1.4, 
U = 48. 

In Fig. H the layer magnetizations m a as well as the 
mean magnetization m = ^ ^ Q m a are shown as a func- 
tion of temperature T. While symmetry requires the dou- 
ble layer film to be uniformly magnetized, the magnetiza- 
tion shows a strong layer dependence for d > 3. The mag- 
netization curves of the inner layers (and for the double 
layer) show the usual Brillouin-type behavior. The trend 
of the surface layer magnetization (d > 3), however, is 
rather different. Note that mi depends almost linearly on 
temperature in the range T/Tq — 0.7 — 0.9 and for thick- 
nesses d > 4. Compared to the inner layers, the surface 
magnetization decreases significantly faster as a function 
of temperature, tending to a reduced Curie temperature. 
However, due to the coupling between surface and inner 
layers which is induced by the electron-hopping, this ef- 
fect is delayed and a unique Curie-temperature for the 
whole film is found. The Curie temperature Tq increases 
as a function of the film thickness d and saturates already 



for film-thicknesses around d — 3 — 5 to the corresponding 
bulk value. AjSimilar behavior was found for a bcc(llO) 
film geometry^. 

The critical exponent of the magnetization (Fig. ||) is 
found to be equal to the mean-field value (3 = h for all 
thicknesses and all other parameters considered. This 
clearly reveals the mean-field type of our approximation 
which is due to the local approximation for the electronic 
self-energy. 
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FIG. 3. Temperature dependence of the layer magnetiza- 
tion m a for a three-layer film on a logarithmic scale in the 
critical temperature range. The critical exponent of the mag- 
netization is given by f3 = 1/2 for all layers (dashed lines). 
Further parameters: n = 1.4, U = 48. 

The layer-dependent quasiparticle density of states 
(QDOS) is shown in Fig. || for three temperatures 
TjT c = 0.1, 0.9, 1.0. Two kinds of splittings are ob- 
served in the spectrum. Due to the strong Coulomb re- 
pulsion U the spectrum splits into two quasiparticle sub- 
bands ("Hubbard splitting") which are separated by an 
energy of the order U . In the lower subband the elec- 
tron mainly hops over empty sites, whereas in the upper 
subband it hops over sites which are already occupied by 
another electron with opposite spin. The latter process 
requires an interaction energy of the order of U. The 
weights of the subbands scale with the probability of the 
realization of these two situations while the total weight 
of the QDOS of each layer is normalized to 1. Therefore, 
the weights of the lower and upper subbands are roughly 
given by (1 — n Q - CT ) and n a - a respectively. This scaling 
becomes exact in the strong coupling limit (U/t 3> 1). 
Since the total band occupation (n ~ 1.4) is above half- 
filling (n — 1), the chemical potential /j, lies in the upper 
subband while the lower subband is completely filled. 
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FIG. 4. Spin-dependent quasiparticle densities of states for 
the outer and inner layer of a d=3 film. Solid lines: majority 
spin direction (a =1); broken lines: minority spin direction 
(a =|). Further parameters: n = 1.4, U = 48. 

For temperatures below the Curie temperature Tq 
an additional splitting ("exchange-splitting") in majority 
(<7 =|) and minority (a —\) spin direction occurs, leading 
to non-zero magnetization m a = n a ^ — n a i- For T = 
the majority QDOS lies completely below the chemical 
potential, the system is fully polarized (n Q | = 1). Thus 
the low-energy subband of the minority spin direction 
disappears and the minority QDOS is exactly given by 
the BDOS of the non-interacting system (Fig. due to 
vanishing correlation effects in the a —\ channel. The 
reduced surface magnetization at T = (see Fig. ||) is, 
therefore, directly caused by the layer-dependent BDOS. 

We like to stress that the spin-splitting does not de- 
pend on the size of the Coulomb interaction U as long 
as U is chosen from the strong coupling limit. Contrary 
to Hartree-Fock theory the spin splitting saturates as a 
function of U for values of about 2 — 3 times the band- 
width of the non-interacting system. The same holds for 
the Curie temperature Tyi-I. 

The temperature behavior of the QDOS is governed 
by two correlation effects (Fig. |]). As the temperature 
is increased, the spin-splitting between a =| and a =\ 
spectrum decreases. This effect is accompanied by a re- 
distribution of spectral weight between the lower and up- 
per subbands along with a change of the widths of the 
subbands. For T = 0.9 Tc one clearly sees that in the 
minority spectrum weight has been transferred from the 
upper to lower subband which has reappeared due to 



non-saturated magnetization at finite temperatures. In 
the a =| spectrum the opposite behavior is found. For 
all T < Tc the spin splitting is significantly larger in the 
inner layer compared with the surface. At T = Tc the ex- 
change splitting has disappeared whereas the correlation- 
induced Hubbard splitting is still present. 

Let us consider the question why the surface magne- 
tization shows a tendency towards a reduced Curie tem- 
perature as seen in Fig. |[ This effect can be understood 
by the above mentioned moment analysis of the BDOS 
(Fig. Q): From bulk systems it is known that an asym- 
metrically shaped BDOS is favorable for the stability of 
ferromagnetism in the Hubbard model. In particular, the 
Curie temperatures increase with increasing skewness of 
the BDOSEl. The same trend shows up in the present 
film-system where the skewness of the BDOS is higher 
for the inner layers compared to the surface-layer. Note 
that this argument is somewhat more delicate than in the 
case of Heisenberg films where the reduced surface mag- 
netization is directly caused by the reduced number of 
interacting sites. 



0.6 

o 

I 0.4 

CD 

| 0.2 



0.0 



I 1 I 1 I 1 I 1 I 
P.ITL 



m=0.6 



T=0.1 T r 



i 



i 1 i 1 i 1 i 1 i 




T=0.8 T r 



i.i. 



i 1 i 1 i 1 i 1 i 



a m 




T=0.9 T r 







2468 10 246810 246 810 
film thickness d 



FIG. 5. Surface-layer magnetization mi, central layer mag- 
netization m c and mean magnetization m as a function of 
the film thickness d and for three reduced temperatures 
T/T c = 0.1, 0.8, 0.9. Further parameters: n = 1.4, U = 48. 

In Fig. H the difference between surface magnetization 
mi, central-layer magnetization m c , and mean magneti- 
zation m is analyzed in more detail. For all film thick- 
nesses where this distinction is meaningful (d > 3), the 
surface magnetization is reduced with respect to the mean 
magnetization. This holds not only for very thin films 
(d — 3 — 5, see also Fig. |[) where some oscillations are 
present that are caused by the finite film thickness, but 
also extends to the limit d — » oo where the two surfaces 
are well separated and do not interact. The surface and 
central-layer magnetizations already stabilize for thick- 
nesses around d = 6. Further, Fig. || clearly shows that 
the reduction of mj drastically increases for higher tem- 
peratures. For T = 0.9 Tq the surface magnetization is 
reduced to about half the size of the magnetization in the 



center of the film. 

The charge transfer due to differing layer occupations 
n a = n a f + n a \, is found to be smaller than w 0.03 at 
T = and is almost negligible for finite temperatures. 

We now like to focus on the magnetic anisotropy en- 
ergy within the model system (|^). The second-order 
anisotropy constant K 2 = ifa.so + ifa.dip is calculated 
via Eq. ( |i"l|) which needs as an input the temperature- 
dependent layer magnetizations of the Hubbard film. The 
dipole constants $ Q for an fcc(100) film geometry are 
found to be: $ = 0.7624 • 4tt/3, $i = 0.1206 • 4tt/3, 
$ 2 = —0.0020 • 47r/3, and are set to zero for a > 3. 

To simulate both, surface and volume contribution of 
the spin-orbit induced anisotropy, we choose the effective 
anisotropy field (|j) in the surface-layer to be different 
from its value in the volume of the film: 



Pa 



Ps 
Pv 



for a 
else. 



l,d 



(19) 



In the perturbational approach only the ratio [3 a /uJo is 
important. Thus we are left with only two parameters 
(fixed at T = 0) to model the RT. 
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FIG. 6. Second-order anisotropy constants for a three-layer 
film in the Fe-type configuration as a function of the reduced 
temperature T/Tc- Further parameters: n = 1.4, U = 48. 

In principle these constants could be taken from exper- 
iment or from theoretical ground-state calculations. In 
our opinion, however, this would mean to somewhat over- 
judge the underlying rather idealized model system. We 
are mainly interested in the question whether a realistic 
temperature-driven RT is possible at all in the Hubbard 
film and by what mechanism it is induced. Therefore, we 
choose the parameters (3s and [3v conveniently, guided 
however by the experimental findings in the Fe-type and 
Ni-type scenarios described in the Introduction. 

Fig. ^ and Fig. |t] show that for appropriate parame- 
ters [3s, Pv both types of temperature-driven RT can be 



found within a three-layer film. The same is found for 
any film thickness d > 3. 

In the Fe-type situation (Fig. ^) we consider a strong 
positive surface anisotropy field (3 s /^o — 3.57T together 
with 0v = 0. At low temperatures the system is mag- 
netized in out-of-plane direction. As the temperature in- 
creases, K%, S o decreases faster than —K 2 , dip because of 
the strong reduction of the surface anisotropy and the 
magnetization switches to an in-plane position. In prin- 
ciple this kind of RT is possible for all film thicknesses 
d > 3. For thicker films (3 s has to be rescaled propor- 
tional to d to compensate the increasing importance of 
the dipole anisotropy. 
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FIG. 7. Second-order anisotropy constants for a three-layer 
film in the Ni-type configuration as a function of the reduced 
temperature T/Tc- Further parameters: n = 1.4, U = 48. 

The Ni-type RT from in-plane to out-of-plane magne- 
tization for a three layer system is obtained by a positive 
volume anisotropy field (3y = 6-7T and a negative surface 
anisotropy fig = — 0.6-7T. At low temperatures the dipole 
anisotropy as well as the negative surface anisotropy field 
lead to an in-plane magnetization. For higher temper- 
atures, however, where the surface anisotropy becomes 
less important because of the reduced surface magneti- 
zation, the positive volume anisotropy field forces the 
magnetization to switch to an out-of-plane direction. 
The ratio (3s /(3v determines for what thickness d this 
type of temperature-driven RT is possible and scales like 
Ps/Pv ~ 1/d for thicker films. 

Note that for both types of RT the values of (3s and 
(3v are chosen in such a way that the system is close to 
a thickness-driven RT. In both cases the RT is mediated 
by the strong decrease of the surface-layer magnetization 
compared to the inner layers as a function of temperature. 

Finally, we consider the dependence on the band occu- 
pation n. In Fig. ||the ratio m\/m c of a three-layer film 
is plotted as a function of the reduced temperature T/Tc 



for different band occupations 1.4 < n < 1.7. Above 
n = 1.78 and below n = 1.0 the fcc(100) Hubbard film 
does not have ferromagnetic solutions whereas between 
n = 1.0 and n — 1.3 there is a tendency towards first or- 
der phase transition as a function of temperaturetH and 
no realistic RT is possible. 
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FIG. 8. Ratio mi/m c between surface- and center-layer 
magnetization as a function of the reduced temperature T/Tc 
for a three-layer film and different band occupations n. In ad- 
dition mi /m c calculated within Hartree-Fock theory is shown 
for n = 1.4. Further parameter: U — 48. 

The ratio mi/m c (Fig. ||) decreases as a function of 
temperature for all band occupations. For n = 1.4 
which has been considered in Fig. || and Fig. [tJ it 
changes from m\/m c = 0.9 at T — to mi/m c = 0.4 
close to Tc- However, for higher band occupations this 
strong temperature-dependent change in m\/m c dimin- 
ishes. From the discussion above it is clear that this 
is unfavorable for the RT. We can thus conclude from 
Fig. H that the possibility of a temperature-driven RT 
sensitively depends on the band occupation n. In ad- 
dition we plotted in Fig. ^| the ratio m\/m c calculated 
within Hartree-Fock theory. The result is completely 
different. At low temperatures we find m\/m c < 1 
whereas mi/m c > 1 close to Tc- Note that the ratio 
mi/m c is almost constant for the wide temperature range 
0.1 < T/Tc < 0.9. Therefore, a realistic temperature- 
driven RT is excluded within the Hartree-Fock approx- 
imation. This holds for all parameters that have been 
considered here. 



VI. CONCLUSION 

We have applied a generalization of the spectral- 
density approach (SDA) to thin Hubbard films. The 



SDA which reproduces the exact results of the t/U- 
perturbation theory in the strong coupling limit, leads to 
rather convincing results concerning the magnetic prop- 
erties. The magnetic behavior of the itinerant-electron 
film can be microscopically understood by means of the 
temperature-dependent electronic structure. 

For an fcc(lOO) film geometry the layer-dependent 
magnetizations have been discussed as a function of tem- 
perature as well as film thickness. The magnetization 
in the surface layer is found to be reduced with respect 
to the inner layers for all thicknesses and temperatures. 
By analyzing the layer-dependent QDOS this reduction 
can be explained by the fact that in the free BDOS both 
variance and skewness are diminished in the surface-layer 
compared to the inner layers. 

The inclusion of the dipole interaction and an ef- 
fective layer-dependent anisotropy field allows to study 
the temperature-driven RT. The second-order anisotropy 
constants have been calculated within a perturbational 
approach. For appropriate strengths of the surface and 
volume anisotropy fields both types of RT, from out-of- 
plane to in-planc (Fc-type) and from in-plane to out-of- 
plane (Ni-type) magnetization are found. For the Ni-type 
scenario the inclusion of a positive volume anisotropy is 
necessary. The RT in our itinerant model system is medi- 
ated by a strong reduction of the surface magnetization 
with respect to the inner layers as a function of tem- 
perature. Here a close similarity to the model calcula- 
tions within Hciscnberg-type systems is apparent, despite 
the fact that these models completely ignore the itiner- 
ant nature of the magnetic moments in the underlying 
transition-metal samples. Contrary to Heisenberg films, 
the band occupation n enters as an additional parame- 
ter within an itinerant-electron model. We find that the 
possibility of a RT sensitively depends on the band occu- 
pation. The fact that no realistic RT is possible within 
Hartree-Fock theory clearly points out the importance of 
a reasonable treatment of electron correlation effects. 
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